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ABSTRACT 

In this thesis the generation of equivalent networks by means of a 
congruent transformation using a variable transformation matrix applied 
to the admittance matrix of an n+l node, common datum, network is 
studied. The ranges of the values of the elements of the transfor- 
mation matrix which yield passive RLC components for the equivalent 
network without mutual inductances are considered. When the trans- 
formation matrix is diagonal the solution involves a set of linear 
algebraic inequalities. In the general case, when all the elements 
of the transformation matrix are present, the solution involves non- 
linear algebraic inequalities. The region of passive component 
realizability consists of a hyper-volume in the Space of the elements 
of the transformation matrix, and the use of exponential searching 
functions to generate equivalent networks with design constraints 


is investigated. 
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Equivalent Networks - Historical Background 


Pe introduction 

It is well known that there is no unique solution to a synthesis prob- 
lem for realizing a given driving point or transfer function. A group of 
networks can always be constructed which have the same input immittance, 
or transfer function, or both. These groups of networks are called 
equivalent networks for the particular network function of interest. In 
general, equivalent networks can be generated by two distinct approaches. 

(a) Using different synthesis procedures for realizing the network 
function, such as Brune, Darlington, Bott-Duffin, Miyata or other pro- 
cedures for the case of input impedance synthesis. These procedures 
are well covered in the literature for input as well as transfer function 
synthesis, Ginllemin [3], Van Valkenburg [8]. 

(b) Using a proper linear transformation of the admittance or im- 
pedance matrix of a given network such that the immittance between any 
desired node pair, or the transfer function relating the response in any 
node pair to the excitation at another pair, is kept invariant. 

The second method is a more general approach and is readily applic- 
able to computer design techniques. It leads to a technique for gener- 
ating a group of equivalent networks among which the best suited to the 
problem at hand can be chosen. 

The method of generating equivalent networks by linear transfor- 
mation was first introduced by Cauer [1] in 1929 in Germany and the 
first English paper on the subject was printed in 1931 by Howitt |2]. 
Guillemin [3], discussing this subject, has shown that the same result 
as obtained by a transformation can be achieved by elementary oper- 
ations on the rows and columns of the network admittance or impedance 
matrix. Schoeffler [4] has shown that the transformation can also be 
applied to the branch matrices, and has also developed a technique of 
generating continuously equivalent networks [5]. 

In this chapter a discussion of the Cauer's, Guillemin and 


Schoeffler's method is presented. 


1.2 Equivalent Driving Point and Transfer Impedance 


The equilibrium conditions of a network containing n+ 1 nodes 
with n ports defined by the terminals between each of the n nodes and 
datum, can be written in the following matrix form: 

Liesl y jiv (aes) 
where: 

[ V ] is the column matrix of the node to datum voltages 

LI ] is the column matrix of source currents connected between each 

node and datum 

LY ] is the n x n node to datum admittance matrix of network 

The admittance matrix is formed in the following manner assuming 
no mutual coupling in the network. 

(1) Each diagonal element an is the sum of all the admittances 
connected to node i. 

(2) Each off diagonal element ij i#j is the negative of the 
admittance connected between nodes i andj. 

The driving point impedance of a network excited at the input node 


i only is given by 


pampaiss (122) 
i A 
where: 
& = the determinent of the admittance matrix [Y ]. 
A.. = the minor obtained from Aby deleting the ‘i row and 


il 
the ee column. 
Equation (1-2) can be obtained by expanding (1-1) and using Crammer's 
rule. 


It has been shown by Howitt [2 ] that if the admittance matrix 


LY ] undergoes a real non-singular transformation I 


' ip 
Wie ahirihy ih (1-3) 
where [I] is a square real non-singular matrix of the same order as 
[Y ] and T denotes transpose, then | Y'] is the admittance matrix of a 


new network with the same number of nodes as the original network. 


In general the input impedance at any node i and the transfer im- 
pedance between any two ports i and j of the new network is different from 
that of the original network. However Howitt [ 2 ] and Schoeffler [4 ] 
have shown that if the fa column of [I] consists of zeros except for 
the ‘as term which is unity, then the input impedance at port i remains 
invariant under Panceemiation. If two such columns are chosen, the 
driving point impedance at the ports corresponding to the columns, as 
well as the transfer impedance between them, is invariant to the trans- 
formation. Example 1-1 illustrates the invariance of the input impedance 
under transformation (1-3) for a simple 3 node circuit. 

EXAMPLE 1-1 


Consider the resistance network of Fig. l-la. 
] 


: 1 ] 
in 
Fig. l-la 
2 -] 
Ly J=(G ]J= 
-] Z 
oe 2 
in A 3 
Using the following transformation matrix: 
(le = : = where a and b are arbitrary numbers yields 
0 b 
Z 
2(l-a-2a ) -b+2ab 
[G'] = 
-b+2ab 2b° 





* 
An original proof based on matrix partioning is given in appendix I. 


resulting in Fig. 1-lb where 


b(1-2a) 
Za ah ae Cimced 
in? 
2 (1 -aua- eee 
Fig. l-lb 
Ae 2b- ; 
Z = —— = 2 2 2 = 
in A 2b (l-ata’') - (-b+2ab) 3 
aw ll 
(a) Substituting a = 3, b=- 7 will lead to the circuit of 
4755 
Zi, a3 4/9 8/11 
in 
rig. lke 


(b) Substituting a = 3, b = -1 will lead to the circuit of Fig. 1-ld. 
WAS 


2 73 


D3 





mig. 1.ld 


As can be seen the 3 networks are equivalent, in that their input im- 
pedance is the same. 

As can be seen from example (1-1), the choice of the elements of 
[I], other than the first column are arbitrary. If any of the elements 
in the transformed network are zero, the number of elements in the new 
network is different from the original network. 

If the network under consideration contains inductance and capaci- 


tances as well as resistances then we can write 


1 -1 
DoS GC eee ae ieee (1-4) 
where 
| Gal = the network conductance matrix 
{C= = the network capacitance matrix 
ee = the network inverse inductance matrix 


By the distribution law of matrix multiplication [ 7 ] the transformation 
[I. ] can be applied to each of the above matrices to obtain the ad- 
mittance matrix of the new network. 


i 


: root] 4+ stc'] (1-5) 


Cy']) = [G') + 


Comparing (1-5) to (1-4) after.the transformation (1-3) 





[G'] = [Alc Wee (16a) 

Cora Limiigen) | 4 (1-6b) 

[is J=0P ll tee (1-6c) 
EXAMPLE 1-2 1 
l 

We } 
Fig. 1-2a 
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Consider the R, L, C network of Fig. 1-2a, with the input impedance 


Secs 


3 449° +3841 


Z(s) = 
35 


Keele 7 
CE = 


oe é 0 


0 i 

Takeses i) = : é 
eee 

Then: 

fate ee 


i 
—_ ae a 
I 
Nm Ww 


-l 


CL Ee Morar 


J 


iH] 
| 
| 
tho fo 


[C'] bor) Pailin 


i 
| 


The new network is shown in Fig. 1-2b. 


changed. 


: 1/2 


Figs 


10 


l=-2b 


The in put impedance is un- 





It can be seen from example 1-2 that the new network has one more 
element than the original. The circuit configuration and type of elements 
connected between corresponding nodes is also different. 

Guillemin [3 ] has shown that an equivalent network can be pro- 
duced by multiplying the same rows and columns of the admittance matrix 
by the same constant. He has also shown that by uSing a technique 
which seems to be similar to linear programming, an equivalent network 
having all positive elements may be produced from a network with a 


negative element. The following example will illustrate the above 


points. 
EXAMPLE 1-3 
Consider the resistive network of Fig. 1-3a. 
(1) 20 ai ) 
1/4 -1/5 1/10 
Fig. 1-3a 


The network admittance matrix is 


60 =210 0 
Le JS -20 25 = le 
0 =16 20 





the input impedance Z. = ohms. 


40 
If we multiply the second row and second column by 6/5 the new ad- 


mittance matrix is: 


60 =A 0 
ae.) = =24 36 =i 
0 =i 20 


dl 


The new circuit is shown in Fig. 1-3b. 


1 t/24 0 PAL 2 3) 
1/36 1/8 
Fig. 1-3b 
= aa ohms 


A: 40 

The following example will illustrate the invariance of transfer 
impedance. 
EXAMPLE 1-4 


Consider the network of Fig. 1-4a, connected to a current source. 
] l 


SY O90O009Q0O O8 


Gy a l 1 1 | 


Fig. 1-4a 


The resistances are in ohms, capacitances in farads and inductances 
in “*henries. 
0 0 0 0 0 0 2 -l 0 
LG j=|0 2 =1 [C]=l0 1 0 Ot \=|=1e- el. 0 
0 -l 2 0 0 ] 0 0 0 


The transformation matrix [I ]is arbitrarily chosen as: 


1 ¢ 0 
[rj=|o0 3 0 
4 
| 0 0 l 
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Then: 


LG']+0 LCc']=0 


NIW NO o 
NN NIW oOo 
oS AI|o © 
oO 
om | 
t 
CO BIO Ww 


The circuit diagram of the new network is shown in Fig. 1-4b. The 


transfer impedance has remained invariant. 


DYES By 3 
— i 


(t)l & 2 4/3¢ 


id |co 


Ws l 


Fig. 1-4b 


1.3 Realizability 
It has been shown by Guillemin [3 ]that if the transformation 


matrix [I] is real and non-singular, the positive definite character 
of the stored energy functions and loss function of the network being 
transformed remains unchanged... Consequently, under the condition 
that [T°] is real and non-singular, a passive network under the trans- 
formation of (1-3)is always transformed into another passive network. 
J.D. Schoeffler [5 ] has demonstrated that the transformed network 

can always be realized by imbedding the original network in an array 

of interconnected ideal transformers. However, the realization of 
transformed network as a transformerless network with positive elements 
only presents particular problems and imposes some restrictions on the 
values of the elements of the transformation matrix [I]. These re- 
strictions and allowable ranges for the values of the elements of [I] 
such that the transformed network can be realized with positive R, L, 

C elements and no mutual inductance is discussed later. 


After the transformation is performed and the admittance matrix of 


13 


the new network is obtained, it is possible to determine whether the 
new network can be realized by positive R, L, C elements. The follow- 
ing are the necessary and sufficient conditions for [ Y'] to be realiz-— 
able with positive R, L, C elements, Howitt [2 ]. 

(1) The off diagonal elements of each of the matrices [G'], 


[C'] and Pr | has to be less than or equal to zero, i.e.: 


o's, a) for all i #j (1-7a) 
Ci < 0 for all i # j (1-7b) 
Tee =. 20 for abl. i = j (1-7c) 


(2) The sum of all the rows of each of the matrices [G ] , 


=i 
[LC ] and LL §] has to be non-negative, i.e.: 


nN 

phe ol = 0 (1-8a) 
: 1) 

j=) 

nN 

a > 0 (1-8b) 
, 1) 

j=l 

a <j 

i Ae > 0 (1-8c) 
jel” 


Conditions (1-7) and (1-8) can be checked by a computer after the 
transformation is performed. 
1.4 Transformation of the Branch Matrices 

Since the elements of a network are explicitly presented in the 
branch admittance or impedance matrices, it is sometimes advantageous 
to perform the transformation on the branch matrices. 

The branch current and voltage matrices of a network can be 


written as: 


il el 
a Li Se: 
i Mian 
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where J..and \jrepresent the current through and the voltage across 
branch i respectively. The equilibrium conditions of the network can 
then be written as 


CV, 1 


[21 bi (1-6a) 


or 


(1,1 


CY ICv,3 (1-6b) 
where LZ] is the branch impedance matrix formed as follows: 
Each diagonal element Z. is the total impedance in branchi. 
Each off diagonal element a, is the mutual inductance common to 


branches i and j and 


-] 
(t= te (1-7) 
If the loop current matrix of the network is defined by 
Hy 
Crq= Jt 
I 
m 


where L is the current in loopi then [I ] and Li, J are related by 


T 
(1.) = Cal] ‘(r] (1-8) 
where [@ ] is the network tie set matrix formed as follows. 
+] if I,, and]I[, are in the same direction 
Each elementa.. = bi i 
-] if Th and I, are in the opposite direction 
0 if branch i is not included in loop i 


Similarly if the cut set voltage matrix [V ] of a network is de- 


fined by 


[v]= |y, 


lo 


where V. is a voltage with an arbitrary direction assigned to the branch i 
ofa tree of a network, then 

(vy, J=(6)' Cv] (1-9) 
where | 6 ] is the network cut set matrix formed as follows: 


+] if Vii and V, are in the same direction 


b 
0 if branch i is another branch of the tree 
which includes V. 


Each element Bs =<-] if V : and V. are in the opposite direction 


Schoeffler [ 4 ]has shown that if a matrix [fl] is formed such that 


epee) 1 Tibet (1-10) 
where LI ] is the transformation matrix defined in section 1.2, then 
instead of applying [I] to the network admittance matrix, we can 
apply Ll] to the network branch admittance matrix {Gad i 

The branch admittance matrix El of the equivalent network is 


then given by 
[yj] = (Nltyitnl" (1-11) 


The advantage of this method is that the elements of the net- 
works are readily seen from Lyd and the need for transformers in 
the synthesis of the equivalent network is apparent from presence of 
off diagonal elements in LZ] . The disadvantage lies in the fact that 
LTl] has to be formed from equation (1-10) which means that cut set and 
tie set matrices have to be constructed before [1 ] can be formed. 
1.5 Continuously Equivalent Networks 
J.D. Schoeffler [5 ] has developed a method whereby a set of 
equivalent networks is produced as a function of an independent 
variable x. The elements of each member of the set of equivalent net- 
works are functions of the independent variable x. For any value of this 
variable a network equivalent to the original network is obtained. 
Consider that the original [Y ] matrix is Yo . Define a new matrix 
Y(x) whose elements are some functions of x such that 


M0) =" ¥ 
O 
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Choose [such that: (where [ U ] is the unit matrix) 
Le=> Ue B JAx (1-12) 
now 


[(Pity@ier)}? 


Beiocd 


(Lu 1+(€B ]4x) LY@I](LU J+0B ]4x) ! 


= [y¥(x)]+ (Ly)10B ]) +B IL y()1) ax 


+{BIJICY(«) J([B ieee 


If Ax is small then[Y'(x)]will differ.fromlY(x)]by a small amount. 
Define: 


[AY (x)] = [¥Y'(x)] - LY(x)] (1-13) 


then Equation (1-10) becomes: 


LAY (x)] 
Ax 


= [BICY@)]+(Y@)1[B ]+(B ]LY@)][B ]! ax’ 


as Ax = 0 


Ld Y (x) J 1 
eee ae Sa) cee ea) 

If the elements of B are chosen arbitrarily, Equation (1-11) can be 
integrated to yield the elements oflY(x)las a function of x. Schoeffler 
integrates Equation (1-11) numerically and changes the assumed values 
of the elements of B as the integration proceeds. In the numerical in- 
tegration process the values of the elements of B are chosen so that 
the new matrix Y'(x) is dominant.* Thus the elements of B are actually 


a point function of x as the integration process proceeds. 


For a network composed of R, Land C we can write 


[Y@)] = [G@)]+SlCG)]+ = LL) (1-15) 


* 
i.e. the realizability conditions of Sec. 1.3 are satisfied. 
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substituting Equation (1-12) in (1-11) then expanding and equating 
terms of equal powers of S on both sides we obtain: 


ALG) Ie 


= [B][G@) 1+ [G@) JLB )° (1-1 6a) 


al = [BJ][C@®) 1+ [c@® ][B]* (1-1 6b) 


=i] 
tt} ta jtr toy + Crt ite 17 (1-16c) 


which can be shown to reduce to the form 


aLGed 1 ~ tM ile 1 (1-17a) 


ee = [M ]{ C(x) ] (1-17b) 


emaiant J= (mM tio 3 (1-17¢) 


or 


cee 1 = tM IL Y¥@) J (1-18) 


where [M ] is a matrix whose elements are some combination of the 


elements of [B ]. Matrix | M ] can be found by solving the matrix 


equations: 

[BE dleco) I 1GG®) Jip] = [mM Ife@)] (1-19) 
OF 

[B ICY) )+ CY) 1£B]? = LM JLY(x) 1 (1-20) 
Ieee 

[MJ= (B] + Cy) ]£B 1° [yt J (1-21) 


Substituting Y(0) = a in Equation (1-21) yields matrix [M ] which 
can then be substituted in Equation (1-18). After integration this 
yields the admittance matrix [| Y(x) ] of the equivalent networks as 


function of the independent variable x, 


18 


The Diagonal Variable Transformation and Equivalent 


Networks with Constraints 

2.1 Introduction 

In generating equivalent networks using Equation(1-3)or in gen- 
erating continuously equivalent networks by integrating Equation (1-15, ) 
the designer has no criterion for choosing the elements of the transfor- 
mation matrix [I] to achieve a design specification such as changing 
the value of an element connected between two nodes in a given manner 
or eliminating a given element from the network. The designer is not 
even able to choose the elements of [T] such that the transformed 
network can be realized by positive R, L, C elements with no mutual in- 
ductance. In this chapter a simplified transformation matrix, namely a 
diagonal matrix, is chosen for the transformation. The acceptable range 
of elements of [T ] such that the transformed network can be realized 
by positive R, L, C elements only is found. By choosing the elements 


5% 
of [[T ] to be exponential functions of a real variable x, i.e. e r 


group of equivalent networks are generated and it is shown that by 
choosing the values of X's certain constraints can be imposed on the 
network. However as we shall see in Chapter 4, the diagonal form can't 
always generate equivalent networks. 

2.2 The Diagonal Transform tion and Realizability Conditions 


The transformation matrix is considered to be a diagonal matrix with 


positive elements, i.e. 


%11 0 ; 0 
talaga gloat a 0 Yo 9 0 (2-1) 
0 0 Vee 


Lg 


where 


Me > 0 for 1" =e... 1 (2-2) 


The admittance matrix of the transformed network then becomes: 
veel LE ILY J igi (2-3) 
In general 
[y]= (G@]+stcJ+ein'} (2-4) 
Substituting (2-4) in (2¢3) we obtain 
[YJ=(TILG I(VJ+s(ri(c J(r)+S(TI(e utr} 
(2-5) 


Since in general [| Y']will also be composed of R, Land C then 


we can write 


Lh) ames SLO 4). [ii (2-6) 
Comparing (2-6) and (2-5) we can write: 

eel |G ) iid (2-7a) 

Oe Vill Cuels! k | (2-7b) 

[1 tS ea ea (2-70) 


A typical term of transformed matrices [G'], | C'] and eae 


of the form 


g 

By assumption of Equation(2-2)all the y 's are positive. Hence 
the sign of the elements of the transformed matrices [G' ], [ C'Jand 
lpi remain unchanged. If the original network was made up of 
passive elemamls the off diagonal elements of the [G ], [C Jand 
eT matrices were negative. Hence the off diagonal elements of 
the transformed matrices are also negative and one of the realizability 
conditions stated in section 1.4 namely Equation(1~6)is always sat- 
isfied and the circuit elements connected between any two nodes other 
than node to datum, remain positive after the transformation. 


The sum of the ae row of the transformed matrices [G'], [C'] 


and Cras are given by the following equations: 


20 


n 1a) 

S' (G) - aoe = = %i%y (2-9a) 
j=l ie) 
n n 

Sch = Dec = Siva we G 2-9b 

(C) a Ps, % i %j Sa ( ) 

nN nN 

ee ae oe = 

oe ee Sa ii (2-9c) 


It is clear from the method of the construction of the admittance 
matrix that the above sum of the rows corresponds to the value of the 
om connected between node i and datum, i.e. Tig! © os and 
ne 10 


If the second of the realizability conditions that is Equation (1-7) 


is to be satisfied we must have: 


n 

S'(G) = z Ya %y Sy > 0 forall ii (20a) 
j= 
n 

S'(C) = L a %s C,, ae of for all i (2-105) 


n 
St(L ) = . Ya %; Lig > @ for alli (2-10c) 
Equations(2-10 represent a set of constraints on the elements of the 
transformation matrix. Since from Equation. (2-2) y, is always positive 

it can be cancelled from the set of Equations (2-10). Thus in general 
there are 3n linear equations that have to be satisfied, n equations for 
each of the |G' ], [C' ] and cet os matrices comprising [Y']. If 
Equations (2-10) are satisfied with the equal sign then the corresponding 
element connected between the corresponding node and datum is zero. 


For example, if S' (G) is zero then the conductance connected 


= © 
10 
between node 1 and datum is eliminated from the circuit. From the 


Zl 


above we can conclude that there are no set of y's which would satisfy 
all of the 3n equations of (2-10) with the equality sign, since this 

will result in an equivalent network with no node to datum elements, . 
i.e. with infinite input admittance at any port. 

The 3n inequalities of the form (2-10) together with the n in- 
equalities of the form (2-2) are a set of 4n inequalities which must be 
satisfied by any set of y's to assure the fact that the transformed net- 
work is realizable with passive elements only. Obviously these in- 
equalities are all satisfied if all of the y's are unity, since the original 
networks results in this case. 

From the discussion on Chapter 1 it is obvious that to keep the in- 
put imittance between node k and datum constant %, must be equal to 
unity and to keep the transfer impedance between parts k and £ con- 
stant we must choose %, and % to be unity. There are then n-1 or n-2 
remaining y's to be chosen satisfying Equations (2-2) and (2-10). For 
example in the case of transfer function equivalence there are a total 
of (n-2) + 3n = 4n - 2 inequalities to be satisfied by n-2 unknown 
y's. These inequalities are considered in the next section to determine 
if a solution is possible, and if soto ascertain specific numerical ranges 
for the y's over which the inequalities are satisfied. These ranges 
define the domain of the set of y's for which the transformation of (2-1) 
generates networks with passive components for which the transfer 
impedance between ports i and j remains invariant. Obviously a 
solution cannot always be expected since the transformed network is 
confined to the class obtainable by the diagonal transformation matrix 
(2-1). The general case where a complete n x n transformation matrix 
is used will be discussed in Chapter 3. 

2.3 The Domain of the Solutions 

Each of 3n inequalities of Equation (2-10) represents a linear 
constraint on (n-1) or (n-2) y's depending on whether the input im- 
pedance or transfer impedance is to remain invariant under the trans- 


formation. When any of the inequalities (2-10) are satisfied with the 


ZZ 


equality sign the resulting equation represents a hyperplane in (n-2) 
dimensional y - space.* This hyperplane divides the (n-2) dimensional 

y - space into two distinct regions. The points on one side of the 
hyperplane satisfy the corresponding inequality. For example, if the 
inequality T'i6 > 0 is used the equation of the hyperplane is T's. = Q. 
Points on one side of the plane yield o's greater than zero and points on 
the other side yield os less than zero,When all of the inequalities 
(2-10) are considered, ‘the domain of the allowable set of y's is de- 


termined by a hypervolume in »space bounded by a set of at most 


(4n-2) hyperplanes, 3n of which are of the form: 


n 
Yee %, 9 = 0 (2-11) 


and the remaining number in the form 


amin 0 (2-12) 


where i#k andi # £ since ae A a 1 due to the required in- 
variance of the transfer function. Points on the planes representing 
Equation(2-12) are not acceptable since these points violate the required 
non-singularity of the transformation matrix (2-1). 

Any point on the hypervolume described above represents a set of 
y's which would satisfy all of the inequalities of (2-2) and (2-10). 

The realiZability volume may be determined in a systematic manner 
as follows. Solve all the possible combinations of the (n-2) simultaneous 
equations, taken from the set of (4n-2) equations of the form (2-11) and 
(2-12) for a set of (n-2) y's, and check whether these values also 
satisfy the remaining (4n-2) - (n-2) = 3n inequalities which are not 
used in choosing the (n-2) simultaneous equations. Any set ofy's 
which satisfy these conditions is a vertex of the realizability volume. 


The hypervolume is bounded on the sides by any of the hyper surfaces 


* 
For convenience of writing, the case of transfer impedance equi- 


valence is considered. The case of input impedance equivalence follows 
in the same manner. 
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which passes through (n-2) vertices. 
In general a total of 
(4n - 2) ! _ (4n - 2)! 
[ (4n - 2) = (n - 2)]! (n - 2) ! 3n! Green): 


linear simultaneous equations must be solved and checked in the above 


(2-13) 


manner. In practice many of the solutions are trivial and the process 
may be programmed on a digital computer. After all vertices of the 
realizability hyper-volume are determined, the boundary hyperplanes may 
be determined from the simultaneous equations whose intersections 
determined the vertices. The point where all of the v's are equal to 

one always satisfies the realizability conditions and is always in- 
cluded in the realizability hyper-volume. 

If there is no set of y's which would satisfy the Equations (2-2) 
and (2-10) we can deduce that the transformation of (2-1) will not pro- 
duce any equivalent network with passive components, and the real- 
izability hyper-volume in this case is reduced to a single point where 
all of the y's are unity. If a realizability volume does exist it is possible 
to remove at least one node to datum element by choosing a set of y's 
represented by a point on one of the boundary hyperplanes of the form 
(2-11). 

EXAMPLE 2.1 
Consider the RLC network shown in Fig. 2.1, 


l ‘ 2 


70m a vA 





Pig .s2am 


where all the conductances are in mhas, all the capacitances in farads 


and all the inductances in henries, 
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The transformed matrices [G'], LC']Jand ee aac 


2 “Yo 9 0 
7 2 
%2 99 0 
LG']= 
: : 33 
0 0 0 
2 
0 
[o'l= |, 
0 


0 
0 
0 
1 
0 0 0 
%92 “¥407%33 
2 
%9%33 9733 £ 
0 0 2 


NS 


2 


rly. ' %22 a 
2 
: %2733 4733 36 
0 0 “Y39 l 


The inequalities (2-10) in this case are 12 inequalities but only 8 of 
them contain the variables Yo and V3 and of these only 6 are distinct and 
are tabulated below. 

From [| G'] we have 


(a) 2 - > 0 


Yo 
(b) -1 + Yop 0 
From [| C']we have 
(C) Yq - ¥33 


> 0 


From aa we have 
(e) “Yo9 + 4Y3 3 = 
- = 
(f) “yz, +1 2 0 
In the y-plane the above inequalities define regions as shown in Fig. 2.2. 
Taking into account the Equation(2-2)which states 


Yoo 7 9 = 6) 


uae 
the realizability hypervolume in this two dimensional example will be a 
four sided polygon shown shaded in Fig. 2.2. 

To solve the above problem with a digital computer we must solve 
all the combinations of the six equations obtained from the inequalities 
a, b, c, d, e and f two ata time. 

Table 2.1 presents all the possible solutions of the group of two 
simultaneous equations together with the letter representing the solutions 
on the y-plane of Fig. 2.2. The solutions which yield valid inter- 
sections which satisfy all of the other inequalities are indicated by the 
asterisks and are the vertices of the realizability area. The boundaries 


of the volume are given by the lines representing a, b, e and f. 
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TABLE 2.1 


EQUATIONS SOLUTIONS INTERSECTION 
PAIR Yo 9 Y3 3 POINT 
ab none none none 
are 2 2 A 
ad Z 2 3 B 
a ex 2 3/4 C* 
ene 2 ] De 
bre* ] 1 E* 
bd ] LZ3 Ie 
b e* l / 2 G* 
b {* ] l E* 
cd 0 0 0 
ce 1/3 1/3 H 
ot * ] ] le 
de 3 ] 

a f . ] 
ert 3 } 
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2.4 The Exponential Searching Function 


If a volume of realizability does exist then for each point in this 
volume there exists an equivalent network. These equivalent networks 
can conveniently be investigated by letting the elements of our trans- 
formation matrix |I ] to be exponential functions of an independent 


real variable x. Thus let 


¥. = ei (2-14) 


where A is a constant which may be chosen arbitrary or in accordance 
with some additional constraint placed upon the desired network. 

The exponential function of Equation (2-14) satisfies the require- 
ment of Equation (2-2) in that it is greater than zero for all finite 
values of the independent real variable x. Hence the transformation 
matrix [I] remains real and non-singular for all finite values of x. 
Moreover when x= 0, %4 = ], and for this value of the variable x the 
original network is reproduced. 

When the exponential functions of Equation (2-14) are used as 


elements of the transformation matrix, [I ] becomes 


ai 0 s : ; 0 
0 BE ; 0 
inc) | = ; ; ; ; QO | (2-15) 
Ann* 
} 0 0 F : , e j 


A typical element connected between two nodes in the transformed 


network becomes: 


(A, + )..)x 
oa = g e JJ (2-16) 
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And a typical node to datum element becomes: 


n n ( ha sme lp 

ge = l=, ore e J (2-17) 
j=0 jel” 

And a typical equation of the hyperplane used in determining the real- 


izability volume becomes: 


y MS 
OQ 
@)) 
1 
© 


(2-18) 


g ij and Tso can easily be plotted by a digital computer as a function 
of x. The range of x for which all of the 3n equations of the form of 
Equation (2-17) remain non Zero is the range for which transformation 
of (2-15) yields equivalent networks with positive R, L, C elements. 


Differentiating (2-15) with respect to x yields: 


otPeol 2 trital= (altri (2-19) 
where: 
Av4 : 
0 doo 0 
[aA] = | ye: (2-20) 
0 0 D 


nn 


Differentiating (2-3) with respect to x yields: 


dt y'] cla 
dx dx 


_ str) 
= (v¥IIP]+(Fily] (2-21) 
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Substituting (2-19) in (2-21) we obtain: 


UY = tity) + tvital (2-22) 
Thus the use of the exponential search function for generating equivalent 
networks is the same as integrating Equation (2-22) with the initial 
condition that LY']=[Y ] when x=0. 

Comparing Equation (2-22) with Equation (1-11) we observe that 
the two are identical with [B ] replaced by [A]. 

The choice of the X's may be made to impose additional constraints 
upon the transformed network as follows. 

(a) To keep the transfer impedance between port k and the output 

voltage between nodes £ and m constant (no common node 


= r = OF 
LL nn 
(b) To keep the admittance between nodes i and j invariant under 


between input and output) set Ar = 


the transformation set dei = “AS . From Equation (2-16) it 
can be seen that in this case 4 = ore ? 


(c) From Equation (2-16) 
( dr. a ok dA. + x 
tip Sj 


g ij = a 2 (2-23) 


and 


Chea) x 
= g e kk Lh (2-24) 


Dividing (2-23) to (2-24) we obtain 


rae: | aaltme lanadidead 
ij kd ij k 


Hence to keep the ratio of the admittance connected between 
nodes i and j to the admittance connected between nodes k 


and £ constant set CA, + 5) = ( Avy + d og) } 
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(d) Since the circuit elements between nodes i and j increase or 
decrease monotonically with x depending upon the sum 
( a + ; ), the X's can be chosen so that a specified com- 
ponent changes more rapidly or slowly compared with others 
over the range of x, 
The use of the exponential search function and constraint (b) for choosing 


the A's is demonstrated in example 2.2. 


EXAMPLE 2,2 
Consider the circuit of example 2.1. 1 
(1) (2) (3) ee.) 
J Z 
co m) l 1/2 e 1 
Figseees 


It is desired to keep the transfer impedance as well as the admittance 
connected between nodes 2 and 3 constant. 


From Fig. 2.2 it can be seen that since x = 0 corresponds to 


point E, to stay in the area of realizability we have to choose doo and 
33 such that Yoo increases and Y 39 decreases with x. Hence 33 
must be negative anda 22 must be positive, i.e. X59 a “33 = constant 


the value of the constant depends on the desired rate of change of the 
circuit elements. Since at the moment we are not concerned with rate 


of change of elements we choose 


Lge ee 


The | Ge gerC *pane or mae matrices would then become: 
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ae. |= 


me. |= 


[L']= 


— 
a) 


The new circuit elements become: 


J 10 


e 


= ae 


x 


erie 


~ 20 


© 23 


ih 


“40 


33 


ZK Ze 
-2xX = 
; _ UR 
i 20 i/(e - 1) 
Se) 
i 5 a? 
i 39 = 1/(4e = 1) 
| _ Xx 
L 34 7 e 
ee =e / (se ~) 
40 


The variation of the conductances, capacitances and the inverse in- 
ductances with the independent variable x are shown in Figs. 2.4, 2.9 


and 2.6 respectively. 
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Fig. 2.4 Variation of Conductances with x, Example 2.1 
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Variation of Capacitances with x, Example 2.1 
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Fig. 2.6 Variation of Inverse Inductances with x, Example 2.1 
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The General Transformation 

3.1 Introduction 

In Chapter 2 a diagonal transformation matrix was used and it was 
shown that the realizability hyper-volume, when it exists, can be 
determined by the manner described in that chapter. The diagonal trans- 
formation matrix however is limited in application. For example itis 
obvious from Equation (2-8) that if no conductance existed between nodes 
i and j in the original network, none of the equivalent networks gen- 
erated will have a conductance between nodes i and j. But the sim- 
plicity of the transformation (2-3), linearity of the realizability 
equations (2-10), and the ease of applying constraints on the transform- 
ed networks warranted the use of the diagonal transformation. 

In this chapter the full n x n transformation matrix [I] is used. 
The equations for the elements of the transformed networks are obtained. 
The realizability conditions and the allowable range for the elements of 
the transformation matrix [I] so as to yield to equivalent networks 
with positive R, L, C circuit elements is then discussed. 


3.2 The General Transformation and Realizability Conditions 
The transformation matrix IT is chosen to be 


abl 112 %n 
YO] %92 : Yon 

Lidl = (3-1) 
el : : Yan 


A typical element of the conductance, capacitance, and inverse 
iImametancesmatrices |G"), [| C'] and Capea of the equivalent networks 


produced will then be of the form: 
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n 
me 22 
‘J 2,k=l 


Yin” jk Fak oe 


Since the original network is made up of bilateral elements, the 
Dea) pw eeedrand oes matrices are symmetric. Hence Equation (3-2) 


can be written as: 


n n 
—s | , + & ar | 3-3 
a", ey eh = eer © viel oe (3-3) 
k=1 p=l 
q=p+l 
Where i =j, Equations (3-3) yield the diagonal elements of the matrices 
we) tc’ ] and [l'-}. ‘huss 
n 2 n 
fee — I Sia Ei, 5 NRE 7 remem) Z7pemenmen (5) GY cama ie (3-4) 
li Pay kk ‘ik el pq‘ ip’ iq 
fehl 


We recall from Chapter 1 that realizability with passive components 
required that the off diagonal elements of the [G' ], LC' ] and pulang 
be zero or negative and that the sum of terms in any row of the above 
matrices be zero or positive. Thus realizability with passive elements 
in the case of the full transformation matrix [IT ] requires that 

oat < 0 fori # j (3-5) 


nN 
' = ' > = 
cae 2 ae 0 (3-6) 


For ann port network each of the matrices [G' ], [C' ] and 
oo | is a symmetrical n Xn matrix. Hence we have in general 
3n(n-1)/2 inequalities of the form (3-5) and 3n inequalities of the form 
(3-6) to be satisfied. Since the transformation matrix [I] need not be 
symmetrical it contains ae elements. The invariance of the input or 


transfer impedance requires that one or two columns of the transformation 
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matrix [I] be unit vectors (i.e. zeros everywhere except on the diag- 
onal position). Hence if the invariance of the transfer impedance is 
required a total of 3n(n+1)/2 inequalities have to be satisfied by 
enon) = n(n-2) y 's if the transformed network is to be realizable with 
passive components. 

Examination of (3-3), (3-4), (3-5) and (3-6) shows that the in- 
equalities (3-5) and (3-6) which have to be satisfied are non-linear 
functions of the y's. Hence the method used in Chapter 2 for the diag- 
onal transformation matrix is not applicable here, since the equations to 
be solved are not simultaneous linear equations. To be able to apply 
the method of Chapter 2 to the transformation of (3-1) used in this chap- 
ter, we will linearize the inequalities of (3-5) and (3-6) as shown in 
the next section. It must be noted that in the case of transformation 
with the full transformation matrix of (3-1),in contrast to the case of 
diagonal transformation matrix ,the y's can take both positive and neg- 
ative values. 

3.3 Linearization 

We know from the discussion in Chapter 2 that the inequalities 
(3-5) and (3-6) are satisfied when real 1 and in O tor iy. 

(«S This will result in the transformation matrix [I’] to be the 
identity matrix [I ] and the transformation will yield the original net- 
work. 

We can linearize the inequalities (3-5) and (3-6) by expanding the 
matrix [Y'] in a Taylor series expansion pout the point Ya = land 
y.. = 0 i #j as shown below: 


1) 
From (1-3) we have 


[yj = (rityitri? (3-7) 
or from (1-6a) we have 
[@}= (Irj(@](ri* (8-8) 
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Hence 


Wien, = [AlIlc iit SermirG mari! 


nN 
“tacief x ie” [G]{L1+{(Tl{e] 





biden O%% 
n 
x(x lay, (3-9) 
L,j=1 oY Fi 
But 
0 0 
T ; : 
atte A vA = ay; 
oe 
0 0 
Hence 
a ee a 
= or A A 
ae ay, = [a2 Ly =[Ay ] 
j=l 
(3=10) 
AY 1 : : ; : AY ae 
LL 
Since we are expanding [G'] with a Taylor series about the point 
y,, = ma i#~j we can write 
= 1+ : 
Vai | 6. (3-11) 
— + = 
Y 5; 0 aT (3-12) 


4l 


where 6. and i, are incremental change on Yi and Vij , EES 

6 = Ay. and oF - AY; and the transformation of (3-1) with 

6. = a5 = 0 yields the original network. The origin of the &space 
corresponds to [I] being the identity matrix. 


Substituting (3-11) and (3-12) in 3-9 we can write 
n 


Met = Tan ij = x 
1,j=1 


d{T ] 
8y A V5 (3-13) 


Substituting (3-13) in (3-9) we obtain 


[AG] = Jere JIT) + fh) vempuen (3-14) 
Expanding [G'] with the Taylor series around [|G ] yields 
re! ==ra nl a + [AG'] lo (3-15) 


where [G' ] - ancePAG" | I denotes the matrices of | G'Jand 
[ AG'] where ae = 1 and Lat 0 i#j or where oy = 6,, =, 
From (3-8) 
Bees LG | (3-16) 


Since [T ] \ = eT |. = [LU] where | U] denotes the identity 


matrix then from (3-14) we have 


[AG'] |, = [6](G]+ fe vie 7" (3-17) 
Substituting (3-17) and (3-16) in (3-15) yields 


feo] (6h + (6)le) ene a (3-18) 


Hence a typical term of the [G'], [C'] and tae | matrices would 
be 


(3-19) 


Il 
Q 

4 
M5 


(9,55 +95, 8) 


where only the first two terms of the Taylor series are taken. 
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Equation (3-19) will enable us to linearize the inequalities of 
(3-5) and (3-6) and hence apply the method of Chapter 2 to obtain the 
area of realizability and the range of y's for which the transformation 
of (3-1) yields equivalent networks with passive elements. In essence 
this procedure identifies a small region in 6- space, any point of which 
corresponds to a set of y's for which the transformation of (3-1)may be 
applied to yield equivalent networks with passive elements. Once 
this region is defined a searching function such as the exponential 
searching function of Chapter 2 can be used to generate equivalent 
networks whose circuit elements are functions of an independent variable 
Such as X and the complete inequalities of (3-5) and (3-6) investigated 
to yield equivalent networks with passive elements. 

3.4 The Exponential Searching Function 
We use an exponential searching function for the elements of the 


transformation matrix LI ] in the following form: 


ee = (3-20) 


u i #j (3-21) 


i 
— 
1 


Viawe 
1) 


Approximating (3-20) and (3-21) by the first two terms of their Taylor 
series expansion and comparing with Equations (3-11) and (3-12) respec- 
tively we obtain 


aC ie (3-22) 
ii li 


-»..x = 6.. (3-23) 
ij ij 


Equations (3-22) and (3-23) in conjunction with the linearized form 
of Equations (3-5) and (3-6) can be used to choose values of the \'Ss 


and the range of x, as best illustrated by the following example. 
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EXAMPLE 3.1 


Consider the RC network of Fig. 3.1. 
(1) ; (2) 





FAG ee Seal: 


It is desired to find the equivalent networks containing passive elements 
which have the same input impedance as the above network. 


The [G ] and [C ] matrices of the above network are: 


C, 0 1 | 
LC j= [LG l=g 


0 C, -l ] 


Starting with the diagonal transformation matrix 


l 0 
De 


: Y992 


The [G'] and | C'] matrices will be 


II “59 C, 0 
LG']=g , [C'j= , 

a 0 Cry 

Y92 %99 2%22 


The realizability conditions are 
} - Yo9 > 0 

- > 
) + Yoo 2 0 


a > 0 


CoY%o0 = 
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It is obvious that the only value of Yo9 which satisfies the above 
equations is Vo9 = 1] which means that the transformation matrix is the 
identity matrix. Hence the diagonal form of the transformation matrix 
will not generate any equivalent network. 


Now we use the full transformation matrix 


: Pi 
Lit = 
Y99 
The [G'] and [ C'] matrices will be 
1-2 +: ; 7 
Vat? 7 712 Y) Y wae Y 90 
LG'j=g 
; 2 
%92%12 ~ 92 %92 
cl Guen- C 
i Our 2%12%22 
[C']= 
C Cpe 
ee LO, 27.99 


The realizability conditions (3-5) and (3-6) become: 


we 2 
() 9'y9 = 2 ¥yot %1 0+ % 10% 90 ~ % 09) G a 


2 2 


4S 


(c) oF Yop aad < 0 


— 2 
d) Cg =C, + Only o%ont% 2) 2 2 


= 2 
te) Cog O2' Me % aa %22) = | o 


(fF) C14 = CoV 2%202 


The linearized form (3-19) of the above equations are 


1 


ib) 6), +55, 


\V 
on) 


Each of the above inequalities when satisfied with equal signs will be 


presented by a line on 69 56 space as shown in Fig. 3.2. The 


22 
linearized realizability region reduces to the portion of the line 


645 + 6,5 = 0 in the second quadrant of the 6 space shown on Fig. 3.2. 


Using the exponential search function, Equations (3-20)and (3-21) 
become 


4 


| 
eel 
’ 
@)) 


Equations (3-22) and(3-23) become 


22 * 22 


ip 


From the conditions of realizability for this case we have 6 


constant. 


the numerical values of the circuit element. 


=e 


m/e 


The value of the constant determines the rate of change of 


Since we are not concerned 


at this moment about the rate of change of elements we arbitrarily 


choose Avo = 19 =] forx>0Q. 


Hence 


Y92 — 
x 


Yap > le 


Substituting the above values in the expressions for g's and c's yields 


Zio ° 
Sy) 
ees 


1 = — x 
Cho C,+C,(1 e) 

q = x 
C0. we 

ts Re ox 
GC 12 = Cre (l-e ) 


From the above equations it can be seen that the range of x for which the 


realizability conditions are satisfied is given by the two equations 


x 
= > 
Cc, +C,(01 0 ee 
Clemil= ce) < 
2 
Hence 
| a e* <—ee 
or See Ce 1, 
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0 





The equivalent network is shown in Fig. 3.3 below and the variation of 
the network elements with x are shown in Fig. 3.4 for the case of g=1, 


er Ce (e™-1) 





| ie pean 
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Fig. 3.4 Variation of Circuit Elements of Example 3.1 
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Conclusion and Suggestions for Further Research 


4.1 Conclusion 

In this paper the congruent transformation is applied to the admit- 
tance matrix of an n+l node common datum network to generate equivalent 
networks which have the same input or transfer admittance, It is shown 
that a transformation matrix whose elements are exponential functions 
of a real variable x, Such as BS can be used to generate a grouh of 
equivalent networks whose circuit elements are functions of the variable 
x. The range of variation of the elements of the transformation matrix 
required to generate equivalent networks with passive circuit elements 
is found by solving a number of sets of equations called the realizability 
equations. If a diagonal transformation matrix is used these equations 
are linear but when a full transformation matrix is used these equations 
are nonlinear, and have to be linearized. If a diagonal transformation 
matrix whose elements are exponential functions of the real variable x 
is used, then certain constraints can be placed upon the generated 
equivalent networks by choosing certain relations between the constant 


AX 


factors ()'s) in the exponent of the exponential function e 
4.2 Suggestions for Further Research 

It is believed that the subject of equivalent networks has a great 
potential for application in the field of network synthesis and analysis. 
The following are some ideas for further research on this subject. 

(a) It may be possible to find the general form of a transformation 
which would take a given form of a synthesised network to another form, 
For example it may be possible to find a transformation matrix [Mf] which 
transforms the Foster form of a synthesised network to Cauer's form for 
the same input impedance. 

(b) The effect of the presence of mutual inductance in a network 
on the transformation can be investigated with the aim of transforming 
a network with mutual inductance to a network without mutual inductance. 

(c) Finally, from curves of Fig. 2.4, 2.5 and 2.6 it can be seen 


that if the values of the components of the network vary according to 


ol 


the curves shown, then the transfer impedance and input impedance 


a a a = 


remain invariant, i.e. the multiparameter sensitivity of the input or 
transfer impedance is then zero. Hence, multiparameter sensitivity 


can be studied from the point of view of equivalent networks. | 
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Appendix I 
Invariance of Input Impedance 


It was stated in Chapter 1 that if the column i of the transformation 
matrix [i ] is a unit vector (zero everywhere except for a one at the we 
row) then the input impedance at port i, that is Zs remains invariant 
under the transformation. In this appendix a proof of this fact based on 
matrix portioning is presented. 

Consider an n+l] node network with a current source connected be- 


tween each of the n nodes and the datum node as shown in Fig. I-l, 





Fig. I-1 


The equilibrium equation of the network can be written in matrix 


form as 


Mies = CY LV] (I-1) 


with [I], [Y ] and [V ] as defined in Chapter 1. 
Expanding the matrix Equation (I-1) and using Crammer rule, the 


voltage V. between node i and datum is given by 











4 455 4: in 
im AW wis VAG 2 ec agen ett al (I-2) 
where: 
4 = the determinent of the admittance matrix [Y ] 
mat = the minor formed from Aby deleting row i and column j 


If the network is excited only at node i by a current source I, then: 


1 = 0 180) rs teal WS ec i #j 
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Hence We = ai (I-3) 


i A i 
and 
oath 
= =a 1-4 
i A ee) 
The invariance of the input impedance at port i under the congruent 
transformation 


weeteiey Veer (1-5) 


with column i of LI] being a unit vector would be proved if we can 
show that the ratio 4. / & of Equation (I-4) remains invariant under 
transformation (I-5). 

For simplicity of proof we consider the invariance of the input im- 
pedance at node 1. Since by remambering the nod@és we can call any node 
node 1, there is no loss of generality and the proof would be acceptable. 
for the invariance of the input impedance at any port. 

Consider a real non-singular transformation matrix [IT] where 
column 1 is a unit vector 


1 Y12 : , ; us 


n 
0 Yo9 ; ; : Yo 

LI j= (I- 6) 
: Yn2 ° ' : ‘nn 


The admittance matrix [ Y'] of the transformed network is given by 


[yd=(£lfyitri! (I-7) 


where [Y ] and [Y'] are as defined in Chapter 1. 
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Since the determinant of the product of matrices is equal to the 


product of the determinants [I] then: 


(A een (i-8) 
where A' = determinent of [ Y'] 
and |I| =determinent of [I] 
We can partitim the [T ], Gene , CY J], LY'] matrices as shown 
below.* 
| 
1 1 
Yul Yin 1 %12 “In 
—_-— —|j— — ———— << a i vemes) ogee meee eee ee os 
Onn 
| 
| ‘i 7 | 0 
f 
e 
( 
¥ Wels. f ! 
| 
Yyy fos Yin 1 | 0 ; : 0 
oe a i a> ae 
a1 
| x a 
x ie on 
| { 
Yn | fin! 
q 
(I- 9) 
Performing the matrix operation we obtain 
Ab 
LX'jJ=Llal]lx lla] (I-10) 


Hence the admittance matrix of the new network is given by: 


* 


Yi, i#xj = the admittance connected between nodes i and j 


Ye = the sum of all admittances connected to node i 
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ewe eee ee eee eel ee 


Lon txlbod- (I-11) 
[y']= 


The input impedance at node i of the new network can be written 


from Equation (I-4) to be 


» . ul _ det. (Lo [X Ife)" 
: 1 


N A! (I-12) 


Since determinant of product of matrices is equal to the product of the 


determinants [ 7 ] we can write - 





fs 2 2 
det (LallX]fe])=lal*|x |= lel" 4, 13) 
where |a| = determinant of [a]. 
Since the lst column of [I ] is a unit vector 
IT| = Jel (I-14) 
Substituting (I-14) in (I-13) we obtain 
Ti 2 
det (Lallalfo]) = |T| a, (I-15) 
Substituting (I-8) and (I-15) in (I-12) yields 
Z 
. 
] | r | Z A A 
Hence’ 
2 = Z, (I-16) 


Appendix II 
Computer Program for Example 2.1 


ee JOB0184F sARDALAN 


10 


22 


23 


24 


25 


PROGRAM EQTRANF1 


DIMENSIONX(500 ) 9G10(500) G12 (500) G20( 500) »630( 500) 640(500) C100 
1500) »€20(500) 2€23(500) »€30(500) »€40( 500) »P20( 500) »P23(500) sP30( 50 
2)9P34(500 )»P40(500) sITITLE(12) 5GA2(500) »GA3( 500) 


GM2=l. 

GM3=-le 

X(1)=0,. 

Az=0.8 

DO 10 [21,500 
G12(1)=EXPF(GM2*X(1)) 
G10(1)=2e-Gl2(1) 
G20(1)=EXPF(2e*GM2*X(1))-G12( 17) 
G30(1)=EXPF (2¢*GM3*X(1)) 

G40(I)=le 

ClO(1)=2.6 

C2301) =EXPF ((GM2+GM3) ¥X(1)) 

C20(1) =EXPF (2 e¢*GM2*X( 1) )-C23(1) 
C30(1)=3e¢*#EXPF(2e¢*GM3*X(1) )-C23(1) 
C40(1)22¢ 

P23(1) =EXPF ( (GM24+GM3) ¥X(1)) 
P20(1)=EXPF (2e*GM2*X(1))—-P23(T) 
P34(1)=EXPF (GM3*¥X(1)) 
P30(1)=4e#*EXPF(2e*GM3*¥X(1))-P23(1)-P34( 1) 
P4O(1)=1,.-P34(1) 

GA2 (1) =EXPF (GM2*X(1)) 

GA3 (1) =EXPF (GM3*®X(1)) 
XO141)=X01)4+(A/5000) 

CONTINUE 

DO 22 JI=lel2 

ITITLE( I) 212 

ITITLE( 1) =8HARDALAN 

LABLE=4HG10 

CALL DRAW(5000Xe2G100l sO0sesLABLEsI TITLE 20s02020 020009505 e0eLAST ) 
DO 23 J=lel2 

ITITLE( I) =12 

ITITLE( 1) =8HARDALAN 

LABLE=4HG12 

CALL DRAW(500 9X0G1292s0sLABLEs ITITLE 2000900000900 595s0eLAST ) 
DO 24 J2=ls912 

ITITLE( I) =12 

ITITLE(1)28HARDALAN . 

LABLE=4HG20 

CALL DRAW(5002X2G2002 sOsLABLEs ITITLE 2090200000 92095955909LAST ) 
DO 25 12112 

ITITLE( I) =12 

ITITLE(1) =8HARDALAN 
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26 


ay 


28 


Zo 


30 


31 


32 


a3 


34 
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LABLE=4HG30 


CALL DRAW(5000X9G3092 s0sLABLEsITITLE 20902090 20909595 908LAST ) 


DO 26 1#1lsl2 
ITITLE(1)=12 
ITITLE(1) =8HARDALAN 
LABLE=4HG40 


CALL DRAW(500 9X 9G4093s0eLABLEsITITLE 020902090908 595s00L AST ) 


DO 27 T1912 
ITITLE(T)=12 

ITI TLE(1)=8HARDALAN 
LABLE=4HC10 


CALL DRAW(5009X0C1001s0sLABLEs I TITLE 2020002090009 595s0sLAST ) 


DO 28 IT=lsl2 
ITITLECI)=12 
ITITLE(1) =8HARDALAN 
LABLE=4HC23 


CALL DRAW(500 9X0C2392s0sLABLEsITITLE 90909000 20909 595909LAST ) 


DO 29 [=lsl2 
ITITLEC IT) =12 
ITITLE(1)=8HARDALAN 
LABLE=4HC20 


CALL DRAW(500 oXeC20s2s0sLABLEsITITLE 9000909000909 595008LAST ) 


DO 30 [=lsl2 
ITITLE(I)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HC30 


CALL DRAW(5009XsC30e2s0sLABLEsITITLE 2090209020 0095e5908LAST ) 


DO 31 [=lsl2 
ITITLECI)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HC40 


CALL DRAW(500 9X9C4003 s0sLABLEsITITLE 2090009090909 5e590eLAST ) 


DO 32 I=1lsl2 
ITITLE( 1) =12 
ITITLE(1)=8HARDALAN 
LABLE=4HL23 


CALL DRAW( 500 oXeP239l1ls0sLABLEsI[TITLE 9090009090909 595e00LAST ) 


DO 33 I=1lel2 
ITITLE(T)=12 
ITITLE(1) =8HARDALAN 
LABLE=4HL 20 


CALL DRAW( 500 sXsP2092s0sLABLEsTITITLE sO 202020 e0s005e5e0eLAST ) 


DO 34 I=1lsl2 
ITITLE(T)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HL34 


CALL DRAW( 500 2XeP3492 0 esLABLE,ITITLE 0202000020909 595e08LAST ) 


DO 35 I1=lsl1l2 
ITITLE( IT) =12 
ITITLE(1) =8HARDALAN 


Site, 
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LABLE=4HL30 

CALL DRAW(500 eXeP3002 s0eLABLEsITITLE 90909090 20s09595n0eL AST ) 
DO 36 J#lsl2 

ITITLE(])=12 

ITITLE(1)=8HARDALAN 

LABLE=4HL4O 

CALL DRAW(500 XeP40e320eLABLE,ITITLE 2090900000 200595009L AST) 
END 

END 
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Appendix Ill 
Computer Program for Example 3.1 


ee J0B0 184F »sARDALAN 


15 


10 


23 


24 


25 


26 


PROGRAM NONLINER 
DIMENSION X(500) 9G10(500) 2612( 500) G20( 500) »C10(500) »C12(500) »C20( 


1500) sITITLE(12) 


A=le 

GM12Z=le 

GM222le 

G=le 

Cl=le 

C2zle 

X(1)=#0.6 

DO 10 [=1s500 

G12(1) =G#¥EXPF ((GM22+GM12)*X(1)) 

C1201) 2C2*#( EXPF(GM22*¥X(1)) ¥CEXPFIGM12*X(1))—-Lle)) 
C1O(1)32C14+C2*(EXPFIGM22¥*X01))*( 1Le-EXPFI(GM1L2*X(1)))+01e¢-EXPFI(GM12*xX 


1(1)))**2) 


C20(1)=-C12(01)+C2*¥EXPF (2¢e*GM2Z2*X(1)) 

PRINT 159X(01)sG12(1) sC1lO(I)sC12(1) sC20(1) 

FORMAT (1X95F10¢6) 

XC I+1)2X(1)+(A/500e¢) 

CONTINUE 

DO 23 IFlsl2 

ITITLE( I) =12 

ITITLE( 1) =8HARDALAN 

LABLE=4HG12 

CALL DRAW(5009XeG129l] s0sLABLEsITITLE 20909090 202009595 9089LAST ) 
DO 24 I=lsl2 

ITITLE( 1) 212 

ITITLE( 1) =8HARDALAN 

LABLE=4HC10 

CALL DRAW(5009XsC1092s0sLABLEsITITLE 290909202020 909 595 900LAST ) 
DO 25 I=1lsl2 

ITITLECI)=12 

ITITLE(1)=8HARDALAN 

LABLE=4HC12 

CALL DRAW( 500 2XsCl2es2s0sLABLEsITITLE 202038090 90909595909LAST ) 
DO 26 I=1ls12 

ITITLE(1I)=12 

ITITLE( 1) =8HARDALAN 

LABLE=4HC20 

CALL DRAW( 500 9X9C2093 s0eLABLEs,ITITLE 20909090 90909595s0sLAST ) 
END 

END 
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